Measurement of the response of a Ga solar neutrino experiment 
to neutrinos from an 37 Ar source 
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An intense source of 37 Ar was produced by the («, a) reaction on 40 Ca by irradiating 330 kg of calcium oxide 
in the fast neutron breeder reactor at Zarechny, Russia. The 37 Ar was released from the solid target by dissolution 
in acid, collected from this solution, purified, sealed into a small source, and brought to the Baksan Neutrino 
Observatory where it was used to irradiate 13 tonnes of gallium metal in the Russian- American gallium solar 
neutrino experiment SAGE. Ten exposures of the gallium to the source, whose initial strength was 409 ± 2 kCi, 
were carried out during the period April to September 2004. The 71 Ge produced by the reaction 71 Ga(v E , e") 71 Ge 
was extracted, purified, and counted. The measured production rate was 11.0^ J} (stat) ± 0.6 (syst) atoms of 
7I Ge/d, which is 0.79^ ™ of the theoretically calculated production rate. When all neutrino source experiments 
with gallium are considered together, there is an indication the theoretical cross section has been overestimated. 



I. INTRODUCTION 

To verify that their efficiencies are well understood, the two 
gallium solar neutrino experiments have measured their re- 
sponse to reactor-produced neutrino sources of known activ- 
ity. Sources of 51 Cr were used in these experiments, as re- 
ported in [1] for SAGE and in [2] for GALLEX/GNO. Table I 
gives some of the relevant details of these experiments and on 
the experiment with an 37 Ar source that is the subject of this 
article. 

There are several advantages of an 37 Ar source compared 
to a 51 Cr source, which we enumerate here. A major ad- 



vantage is that the desired active isotope must be chemically 
separated from the target following irradiation. Although this 
involves an additional processing step, this separation serves 
to remove almost all impurities that are present in the target. 
Thus an 37 Ar source, in contrast to a 51 Cr source, which is 
made by irradiating Cr metal, without any subsequent purifi- 
cation, can be made practically free of radioactive impurities. 
Other advantages of 37 Ar compared to 51 Cr are that the half- 
life is longer (35 d compared to 27 d), thus giving more time 
to prepare the source and to make measurements, that the neu- 
trino energy is greater (811 keV compared to 747 keV), thus 
giving a higher cross section, that the decay is purely to the 
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TABLE I: Comparison of source experiments with Ga. When two uncertainties are given, the first is statistical and the second is systematic. 
When one uncertainty is given, statistical and systematic uncertainties have been combined in quadrature. The values of R for GALLEX are 
from [3]. 



Item 



GALLEX Crl [2, 3] 



GALLEX Cr2 [2, 3] 



SAGE M Cr [1] 



SAGE J/ Ar 



Source production 

Mass of reactor target (kg) 

Target isotopic purity 

Source activity (kCi) 

Specific activity (kCi/g) 
Gallium exposure 

Gallium mass (tonnes) 

Gallium density (10 21 71 Ga/cm 3 ) 

Measured production rate p ( 71 Ge/d) 

R = p(measured) /^(predicted) 



35.5 
38.6% 50 Cr 
1714^° 
0.048 

30.4 (GaCl 3 :HCl) 
1.946 
11.9 ± 1.1 ±0.7 
1 00+° 11 



35.6 
38.6% 50 Cr 
1868 + g 
0.052 

30.4(GaCl 3 :HCl) 

1.946 
10.7 ±1.2 ±0.7 

0-8ltiii 



0.512 
92.4% 50 Cr 
516.6 ±6.0 
1.01 

13.1 (Ga metal) 

21.001 
14.0 ± 1.5 ±0.8 
0.95 ±0.12 



330 

96.94% 40 Ca (natural Ca) 
409 ± 2 
92.7 

13.1 (Ga metal) 

21.001 
11.0+' " ±0.6 
19* om 



ground state (100% compared to 90%), thus giving a mono- 
energetic neutrino source, and that there are no accompanying 
y rays (except for inner bremsstrahlung [4]), thus requiring 
little shielding and yielding a very compact source. Finally, 
since nearly 97% of Ca is 40 Ca, no isotopic separation is re- 
quired before irradiation. In contrast, to make a 51 Cr source, 
the irradiated isotope 50 Cr must be enriched as its content in 
natural Cr is only 4.3%. 
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FIG. 1: Levels in 71 Ge that can be reached by capture in 71 Ga of a 
neutrino from 37 Ar or 51 Cr sources. 



The use of 37 Ar to measure the response of radiochemical 
solar neutrino detectors was originally proposed by Haxton 
[5], who pointed out that the neutrino energy from 37 Ar decay 
is very near that of the principal neutrino line produced by 7 Be 
electron capture in the Sun at 863 keV. A practical method to 
make an intense 37 Ar source by the (n, a) capture reaction on 
40 Ca at a reactor with a high flux of fast neutrons was given 
by Gavrin et al. [6]. 

Both 37 Ar and 51 Cr can excite only the lowest three energy 
levels in 71 Ge, as shown in Figure 1. The cross section for 
the transition to the ground state of 71 Ge is well determined 
by the decay constant; the strength of the transitions to the 
two excited states, at 175 keV and 500 keV, is much more 
poorly known, but can be estimated from (p, n) scattering [7] 
measurements. As evaluated by Bahcall [8], 95% of the cross 
section for the neutrinos from 37 Ar is to the ground state, the 
state at 500 keV contributes 3.6% to the cross section, and the 
state at 175 keV contributes the remainder, although only an 
approximate upper limit for its strength has been measured. 



TABLE II: Decay modes of 37 Ar and the energy released. The 
atomic-electron binding energies are from [9]. The branching frac- 
tions are based on the ratios L/K = 0.0987 ± 0.003 [10] and 
M/L = 0.104+°;^ [11]. The internal bremsstrahlung fraction and 
average energy are estimated for ls-electron capture using Eq. (2) 
and (3), respectively, in [12], and for 2s-electron and p-electron cap- 
ture are taken from Fig. 1 of [13]. These calculations have been ver- 
ified experimentally for 37 Ar in [13] (spectrum) and [14] (absolute 
intensity). See also [15, 16]. 



Decay 
mode 



Atomic energy 
release (keV) 



Fraction of 
37 Ar decays 



Energy per 
37 Ar decay (keV) 



K capture 2.8224 
L capture 0.2702 
M capture 0.0175 
Is int. brems. 325 (average) 
2s int. brems. 325 (average) 
p int. brems. -10 (average) 

Total 



0.9017 + 0.0024 
0.0890 ± 0.0027 
0093 +0 0006 
0.0005 
0.00007 
0.00007 



2.5450 ± 0.0068 
0.0240 ± 0.0007 
0.0002 
0.16 ±0.02 
0.021 ± 0.002 
0.0007 

2.751 +0.021 



II. DECAY OF 37 Ar 

37 Ar decays to 37 C1 with a half-life of 35.04 + 0.04 d [17], 
as shown in Fig. 2. The decay is solely by electron capture 
and the Q-value is 813.5 keV. 

The decay energy that heats the source is given in Table II. 
The major uncertainty in the source heating is from the frac- 
tion of internal bremsstrahlung (IB) decays and is estimated 
to be ~10%. A Monte Carlo calculation showed that the IB 
y rays not absorbed in the shield outside the source take away 
an inconsequential 0.2 + 0. 1 % of the source energy. We thus 
assume the heat deposited in the source following IB is simply 
the average decay energy. 



3?C1 (stable) 



3 7 Ar (35.04 days) 



813 keV v(9.8%) 
811 keVv (90.2%) 



FIG. 2: 37 Ar decay scheme showing the neutrino energies. 
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III. SOURCE PRODUCTION 

The source was made by irradiating calcium oxide in the 
fast neutron breeder reactor BN-600 at Zarechny, Russia. The 
total fast flux at this reactor is 2.3 x 10 15 neutrons/(cm 2 s), 
of which 1.7 x 10 14 neutrons/(cm 2 s) have energy above the 
2-MeV threshold of the production reaction 40 Ca(n, ff) 37 Ar. 
Nineteen irradiation assemblies, each of which contained 
17.3 kg of CaO (12.36 kg Ca), were placed in the blanket 
zone of the reactor. Irradiation began on 3 1 October 2003 and 
continued until 12 April 2004, the normal reactor operating 
cycle. After a cooling period of a week, the assemblies were 
removed from the reactor and moved to a hot cell, where each 
was opened and from which the capsule containing the CaO 
was removed. These capsules were transported to an extrac- 
tion facility at the Institute of Nuclear Materials, where each 
capsule was cut open in a vacuum system and the CaO dis- 
solved in a nitric acid solution. 

37 Ar was extracted from this acid solution by a He purge, 
purified, and then stored on charcoal at LN2 temperature. The 
purification involved flowing the gas over zeolite at room tem- 
perature, followed by two Ti absorbers, operating at 900- 
950 °C and 400-450 °C. The purified 37 Ar, whose volume was 
~2.5 1 at STP, was then adsorbed on another charcoal trap and 
measurements of gas volume and isotopic composition were 
made. 

As the last steps of source fabrication, the purified Ar was 
transferred to a pre-weighed source holder, which consisted 
of a sealable stainless steel vessel with a volume of ~180 ml. 
Inside this vessel was 40 g of activated charcoal onto which 
the purified 37 Ar was cryopumped. When essentially all the 
37 Ar had been adsorbed, the vessel was closed by compress- 
ing two separate knife-edge seals, one onto a copper gasket 
and another onto a lead gasket. The source holder was then 
weighed to determine the amount of 37 Ar contained within. 
The calculated gas pressure in the source holder was ~17 atm 
at room temperature. To complete the source, the source 
holder was placed within two concentric stainless steel ves- 
sels with a Pb shield between them. These two vessels were 
welded shut and the heat output of the finished source was 
measured with a calorimeter. These procedures were com- 
pleted on 29 April 2004 and the source was immediately flown 
by chartered plane to the Mineralnye Vody airport, close to 
the experimental facility at the Baksan Neutrino Observatory 
in the northern Caucasus mountains. 



IV. USE OF THE SOURCE AT BAKSAN 

The experimental procedures and equipment were basically 
the same as used with the SAGE Cr experiment in 1995. 
These were described in detail in a previous article [1], and 
will only be briefly summarized here. 

A diagram of the experimental area is given in Fig. 3. The 
gallium is contained in seven chemical reactors, designated on 
this figure as numbers 2, 3, 4 and 7-10, all of which are fitted 
with stirrers and ancillary extraction apparatus. Reactors 2-4 
contained 22 tons of gallium and were not used in the exper- 



iment reported here; rather, they were used at the time of the 
Ar experiment for solar neutrino measurements which will be 
reported elsewhere. 




FIG. 3: Plan view of the laboratory showing the ten chemical reac- 
tors, irradiation reactor 6 with the adjacent calorimeter, and the pump 
for transferring Ga between reactors. 

Reactor 6 has the extraction equipment removed and re- 
placed by a Zr re-entrant tube on its axis that extends to the 
reactor center. A specially designed remote handling appara- 
tus is able to grasp the 37 Ar source and place it either at the 
center of reactor 6 or in an adjacent calorimeter. Reactors 7- 
10 all contain about 6.5 tonnes of gallium. To make a source 
exposure the gallium from two of these reactors (either 7 and 
8 or 9 and 10) is extracted so as to remove any Ge that may be 
present and pumped into reactor 6. The source is then placed 
at the center of reactor 6 and the exposure begins. The first 
two exposures were for seven days, followed by exposures of 
approximately two weeks duration. At the end of exposure, 
the source is moved to the calorimeter and the Ga is pumped 
back to the original two reactors where the Ge is extracted. At 
the same time, a clean-up extraction of the Ga from the other 
reactor pair is made and its Ga is transferred to irradiation re- 
actor 6 to begin the next exposure. 

When the source arrived, it was removed from the shipping 
cask, a brief measurement of its y activity was made, and it 
was immediately placed into irradiation reactor 6, which had 
been previously filled with 13.1 tonnes of Ga metal. The first 
irradiation of Ga began at 04:00 on 30 April 2004 (local time, 
GMT+3 h), which we call the reference time. 

Ten extractions were made, designated Ar 1 through ArlO. 
The data for each extraction are given in Table III. As a check 
on the extraction process, a second extraction was made after 
exposure 3, designated Ar 3-2. While the irradiation reactor 
was being emptied and then refilled, the source was moved 
to the calorimeter where its heat output was measured, as de- 
scribed in the next section. 

An outline of the extraction procedure is given in [18]. A 
slight variation of this procedure was used for the Ar exper- 
iment: the usual volume of extraction reagents was divided 
into two parts. The first part was added, mixed, and removed, 
and then this process was repeated with the second part. This 
procedure of making 'two half extractions' results in reduced 
heating of the Ga metal, and thus a higher overall efficiency 
when the two extracts are combined. It was introduced with 
the solar neutrino extractions beginning in 1998 and has been 
used ever since. 

The efficiency of extraction was measured by adding to the 
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TABLE III: Extraction schedule and related parameters. The times of exposure are given in days of year 2004. 



Extraction 


Extraction 


Source 


exposure 


Source activity (kCi) 


Solar neutrino exposure 


Mass Ga 


Carrier 


Extraction efficiency 


name 


date (2004) 


Begin 


End 


Begin 


End 


Begin 


End 


(tonnes) 


mass (fig) 


from Ga 


into GeH 4 


Ar 1 


6 May 


121.17 


127.71 


409 


359 


118.48 


127.78 


13.085 





0.93 


0.59 


Ar Z 


14 May 


1Z0.4Z 


Lij. 1 1 


354 


307 


125.38 


135.78 


1 Q HQ/1 


zij 


0.96 


0.93 


Ar3 


29 May 


136.42 


150.71 


302 


228 


133.51 


150.81 


13.063 


211 


0.93 


0.90 


Ar3-2 


30 May 


136.42 


150.71 


302 


228 


150.91 


151.81 


13.049 


274 


0.93 


0.87 


Ar4 


13 Jun 


151.42 


165.71 


225 


169 


147.47 


165.77 


13.055 


208 


0.97 


0.92 


Ar5 


28 Jun 


166.40 


180.71 


167 


126 


162.47 


180.77 


13.018 


210 


0.98 


0.97 


Ar6 


13 Jul 


181.42 


195.71 


124 


94 


173.57 


195.77 


13.025 


219 


0.98 


0.97 


Ar7 


28 Jul 


196.42 


210.71 


92 


70 


193.49 


210.79 


12.974 


215 


0.98 


0.97 


Ar8 


12 Aug 


211.42 


225.71 


69 


52 


208.48 


225.78 


12.997 


209 


0.98 


0.96 


Ar9 


27 Aug 


226.42 


240.71 


51 


38 


223.47 


240.77 


12.945 


214 


0.98 


0.96 


ArlO 


11 Sep 


241.42 


255.71 


38 


29 


238.38 


255.78 


12.969 


211 


0.98 


0.96 



Ga a known mass of inactive Ge carrier before the start of 
exposure to the neutrino source and measuring the mass of 
extracted Ge. The mass of added carrier is given in Table III, 
as is the efficiency of extraction from the metal and the overall 
efficiency including synthesis into the counting gas GeH 4 . 

The carrier consists of slugs of an alloy of Ge in Ga metal. 
Two slugs were added to the gallium in each extraction reactor 
(numbers 7 and 8 or 9 and 10, depending on which pair of re- 
actors was being used) before transfer to irradiation reactor 6. 
The mass of each slug is known and the average concentra- 
tion of Ge is determined by extraction from a large number of 
slugs. 

Some anomalies occurred in the first extraction Ar 1 . First, 
due to a miscommunication, carrier was not added at the start 
of exposure. The extraction efficiency thus cannot be deter- 
mined in the usual way and it is taken to be the same as the 
extraction efficiency from this reactor pair in the next extrac- 
tion. Second, an error was made during GeH4 synthesis which 
resulted in the loss of some of the sample. Since inactive Ge 
was added to the solution before synthesis the amount lost 
could be determined accurately. 

As noted above, a clean-up extraction was made from each 
new batch of gallium shortly before exposure to the source 
commenced. This second extraction served to remove almost 
all traces of Ge carrier that remained because of the ineffi- 
ciency of the first extraction and also removed any 71 Ge pro- 
duced by solar neutrinos during the -2-week interval since 
the last extraction from this batch of Ga. Since there is no 
leftover carrier, the determination of the extraction efficiency 
is simplified, and since there is no remaining 71 Ge, the analy- 
sis can assume no carryover of 71 Ge from one exposure to the 
next. 



V. MEASUREMENT OF SOURCE ACTIVITY 

Three different methods were used to measure the source 
activity as it was fabricated, a series of eleven measurements 
were made while the source was used at Baksan, and two dif- 
ferent measurement methods were used after the source was 
returned to the fabrication facility. In this section we describe 
these methods and give the results. 



A. Measurements at Zarechny during source fabrication 

In the first method, carried out after argon purification, and 
while the gas was put into the source holder, the total vol- 
ume of gas and its isotopic composition were measured. The 
composition was determined with a mass spectrometer and the 
results are given in Table IV. The gas volume was measured 
by warming the charcoal trap onto which the gas had been 
frozen so that the gas expanded into a calibrated volume, and 
reading the pressure. The gas was then frozen into the source 
holder, and the difference in pressure between before filling 
and after filling implied the volume of gas in the holder was 
2.665 + 0.048 1 at STP. Combining this with the isotopic com- 
position and correcting for decay between the time of volume 
measurement and the reference time gives a source activity of 
409 + 5 kCi at 04:00 on 30 April 2004. The stated uncertainty 
has 68% confidence and includes all known systematics. 

TABLE IV: Gas content of the 37 Ar source 47.5 h prior to the refer- 
ence time in percent by volume. The uncertainty shown is statistical; 
there are additional systematic components whose sum is no more 
than 0.8%. 



H 2 


J7 Ar 


j! *Ar w Ar 


40 Ar 


0.26 ± 0.07 


96.57 ±0.13 


1.87 + 0.06 0.35 +0.03 


0.95 ± 0.03 



In the second method, the source holder was evacuated and 
weighed before filling and then weighed again after filling 
with the extracted gas sample. The difference in mass was 
4.400 + 0.042 g at the time of filling (06:25 on 28 April), from 
which the activity is calculated to be 412 + 3 kCi at the refer- 
ence time. 

In the third method, the heat output of the source was 
measured in a heat flow calorimeter with an inner Pb y-ray 
shield and an outer massive container for thermal stabiliza- 
tion. The heat produced by the source was conducted away 
by a thermopile containing 1500 thermocouples connected in 
series. The calorimeter was calibrated using electrical heaters 
of known power and the thermopile EMF over the range of 
(6-8) W (the expected source power) was found to have the 
constant value 9.81 + 0.02 mV/W. After stabilization of the 
calorimeter with the source the average measured thermopile 
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EMF was 65.9+0.03 mV at 20:00 on 28 April 2004. Applying 
a decay factor of 0.9740 gives a power of 6.54 + 0.04 W at our 
reference time. Using the conversion factor in Table II gives 
the source activity at this time as 401 + 4kCi. The uncertainty 
estimate includes the calibration uncertainty, the errors in the 
calorimeter measuring circuits, and the uncertainties in both 
decay energy and 37 Ar half-life. 



B. Calorimetric measurement at Baksan 

The same calorimeter that was used for the SAGE 5I Cr 
source was again used at Baksan to measure the heat out- 
put from the 37 Ar source. This device has been described 
in [1] and the reader is referred to that publication for de- 
tails. The energy released per 37 Ar decay is, however, less 
than one-tenth that per 51 Cr decay and thus modifications to 
the calorimeter to improve its sensitivity were necessary. Four 
of the eight copper heat conductors between the copper cup 
and the surrounding copper can were removed and voltmeters 
more sensitive than originally provided were used to read the 
thermocouple EMF. 

The calorimeter was calibrated using electrical heaters 
made from Al and Fe and with a mock-up source made from 
identical materials as the real source. The temperature re- 
sponse with the Al heater matched that of the mock-up source 
very well and thus the Al measurements were used for the 
calorimeter calibration. 

The calorimeter was used to measure the source power at 
the end of each gallium exposure to the 37 Ar source. Two 
measurements were made after the final extraction ArlO, thus 
yielding a total of 1 1 measurements, which are plotted in Fig- 
ure 4. 

If a weighted fit is made to this data with a decaying ex- 
ponential whose half-life is fixed at 35.04 d the power at the 
reference time is 6.907 + 0.013 W. x 1 for this fit is 1 1.2 with 
10 degrees of freedom (DOF) (probability = 34%). As a 
check, the same fit was made allowing the decay constant to 
be a free variable, along with the power at the reference time. 
The resultant best fit half-life is 34.80 ± 0.20 d, in agreement 
with the known value. ^- 2 /DOF = 9.8/9 for this fit. 

Using the energy release given in Table II and the conver- 
sion factors 1.6022 x 10~ 16 (W s)/keV and 3.7 x 10 I() decays 
of 37 Ar/(Ci s), the inferred source activity at the reference 
time was 423.5 ± 0.8 kCi. The quoted uncertainty here of 
0.2% is solely from the measurement errors. There are sev- 
eral additional systematic uncertainties that must be included 
in a full uncertainty estimate. These include uncertainties due 
to the differences in thermal properties between the source 
and the calibration heaters (estimated to be 1.5%), in the en- 
ergy release (0.8%, see Table II), in the incomplete absorption 
of the inner bremsstrahlung component of the energy release 
(~ 0.2%), in the 37 Ar half-life (0.04%), and in the capture of 
some of the y rays from the source in the outer parts of the 
calorimeter, thus disturbing the 'cold' junction temperature. 
We assign a total uncertainty for this measurement of +2%, or 
+9 kCi. 

The traditional source of uncertainty in calorimetric mea- 




B 6.4- 



£ 1 . , ■ , ■ , ■ , ■ , ■ \ ^-n >■ 

z 20 40 60 80 100 120 140 
Time since 04:00 on 30 April 2004 (d) 

FIG. 4: Source power measurements at Baksan. The solid curve is a 
weighted fit of the data to a decaying exponential with the half-life of 
37 Ar. In the lower panel the power is normalized to the time shown 
and the 68% confidence band for the fit is indicated by the dashed 
lines. 



surements is the contribution to the source heat from impu- 
rities. For the 37 Ar source, the only significant impurity was 
39 Ar, a pure f3 emitter with a half-life of 269 y. The mea- 
sured 39 Ar concentration at the time of source fabrication was 
0.35% by volume. Since each decay gives 565 keV and the 
volume of the source gas was nearly 2.7 1, the heat from 39 Ar 
decay was only 1.8 mW, 0.03% of the heat from 37 Ar at the 
reference time. The influence of impurities on the calorimetric 
determination of the activity is thus negligible. 

C. Measurement at Zarechny by 37 Ar counting 

The 37 Ar source was returned to the fabrication facility in 
December 2004. The source holder was opened with a spark 
discharge in a vacuum system, the entire gas sample was re- 
moved, and samples of the gas were taken for activity mea- 
surement in proportional counters. At this time the 37 Ar had 
decayed by a factor of ~300. 

Because the specific activity was still very high, it was nec- 
essary to make several volume divisions to reduce the count 
rate to a value that was measurable in a proportional counter. 
These dilutions employed He as a carrier gas since the total 
gas volume and thus the pressure was very low. The fraction 
of gas transferred to the proportional counter was only about 
1 part in 10 8 . To be certain that the gas content of the last 
sampling volume was fully transferred to the counter, spe- 
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cial proportional counters were built which contained a side 
arm filled with charcoal onto which the Ar sample was cryop- 
umped. Including the dead volume in the side arm, the count- 
ing efficiency for the K peak of 37 Ar was 59%. Despite the 
great volume division, the count rate was still very high, sev- 
eral hundred thousand per second. This rate was measured 
by continuously recording the pulses from the counter using a 
waveform analyzer with long time span and then counting the 
number of pulses during a selected time interval. 

Five samples were measured in two proportional counters 
using different methods of volume division. Assuming an 
37 Ar half life of 35.04 d, the weighted average of these mea- 
surements gives a source strength at the reference time of 
405.1 + 2.7 kCi where the uncertainty includes the system- 
atic uncertainties from counting statistics, volume division, 
and counting efficiency. Since the time delay from the ref- 
erence time to the time of these measurements was 287 d, this 
result is rather sensitive to the value of the half-life used in the 
decay correction. The 37 Ar half-life uncertainty in the most 
recent data compilation for this nuclear mass [17] is given as 
+0.04 d, which leads to an additional uncertainty in the source 
strength of ±0.65%. 

D. Measurement at Zarechny by isotopic dilution 



strength of 410 + 5 kCi at our reference time. The stated un- 
certainty includes the error in the half-life. 



E. Summary of source strength measurements 



TABLE V: Summary of all source activity measurements. The stated 
uncertainties include all known systematics. 



Measurement 


Activity (kCi J 'Ar at 


method 


04:00 on 30 April 2004) 


Volume of gas 


409.3 ± 5 


Mass of gas 


412.3 + 3 


Calorimetry at Zarechny 


401.3+4 


Calorimetry at Baksan 


423.5 ± 9 


Proportional counter 


405.1 ±4 


Isotopic dilution 


410.1 ±5 



The results of the six activity measurements are summa- 
rized in Table V. Their weighted average is 409 + 2 kCi. The 
two measurements with the calorimeter are each about 1 .5cr 
above and below this average, but still for all measurements 
X 1 = 7.2, which with 5 DOF has a probability of 21%. 



When the source was opened and gas samples were mea- 
sured in proportional counters, additional samples were taken 
for measurement of the volume concentration of the Ar iso- 
topes with a mass spectrometer. The gas from the source was 
diluted by adding a measured volume of Ar gas from the at- 
mosphere, and the concentration of the Ar isotopes in small 
samples was again measured. From the combination of these 
measurements, using the fact that the volume of the samples 
taken for measurement was a small fraction of the total gas 
volume, we calculate the volume of 37 Ar in the source at the 
time the initial isotopic composition measurement was made 
to be 

V37(before) = , (1) 

fl(after)D(decay) - fl(before) 

where /^(before) and /?(after) are the measured ratios of 40 Ar 
to 37 Ar concentration by volume, before and after the dilu- 
tion, respectively; / 40 = 0.99600 is the fraction of 40 Ar in 
atmospheric Ar; and D(decay) is the 37 Ar decay factor during 
the time between the 'before' and 'after' measurements. The 
added volume of air Ar was V(add) = 27.26 ± 0.49 cm 3 at 
STP. 

This method was applied with several different samples 
taken at various times before and after dilution. As an exam- 
ple, for one such measurement, with the 'before' sample taken 
on 15 February 2005 and the 'after' sample one day later, with 
fl(before) = 2.70 ± 0.02 and fl(after) = 6.49 + 0.10, the cal- 
culated volume of 37 Ar in the source on 15 February 2005 
by Eq. (1) was 7.707 cm 3 at STP, equivalent to 1.281 kCi. 
Extrapolating back to the reference time of 04:00 on 30 April 
2004 implies an initial source strength of 408 kCi. A weighted 
average of all the 'before' and 'after' samples gives a source 



VI. 7 1 Ge COUNTING AND EVENT SELECTION 

The extracted Ge was synthesized into the counting gas 
GeH4, mixed with inactive Xe, and inserted into special ul- 
tralow background proportional counters with a carbon-film 
cathode [19]. The counters were measured in an electron- 
ics system which recorded the full waveform for 800 ns after 
pulse onset. Both event energy and rise-time were obtained 
from the pulse waveform and used to select candidate 71 Ge 
events. The system operated very stably with no need for 
equipment replacement during counting. The counting infor- 
mation is given in Table VI. 

The counters were calibrated with 55 Fe and Cd-Se at the 
start of counting and then approximately every 2 weeks un- 
til counting ended. Extreme stability was observed with the 
peak position rarely differing from one calibration to the next 
by more than 1%, the statistics on the number of calibration 
events. The rise-time characteristics were also highly stable. 

The location of the energy acceptance windows for 71 Ge 
events in the L and K peaks was set from the 55 Fe calibration 
in the same way as for solar neutrino runs. This included our 
standard adjustment factors because of non-linearity in the en- 
ergy scale [18]. 

Exactly the same procedures were used to select candidate 
71 Ge events as we use for solar neutrino runs. As the first step 
in event selection, two time cuts were applied to the data to 
suppress false 71 Ge events produced by Rn. To reduce the ef- 
fect of Rn external to the counters which may enter whenever 
the counting system shield is opened, a cut was made to delete 
the first 3 h of data after the shield is closed and counting re- 
sumes. To reduce the effect of Rn internal to the counters 
which may enter when the counters are filled, a second time 
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TABLE VI: Counting parameters. A is the exponentially weighted live time. The live time and A include all time cuts. 





Counter filling 




Operating 


Counting efficiency after 








Extraction 


GeH 4 


Pressure 


Counter 


voltage 


rise time and energy cuts 


Day counting 


Live time of 




name 


fraction (%) 


(mm Hg) 


name 


(V) 


L peak 


K peak 


began in 2004 


counting (days) 


A 


Arl 


6.5 


695 


YCT3 


1124 


0.346 


0.386 


129.1 


104.0 


0.870 


Ar2 


8.1 


660 


YCT5 


1131 


0.354 


0.382 


136.7 


106.7 


0.843 


Ar3 


8.5 


710 


YCT 1 


1164 


0.347 


0.388 


151.7 


138.0 


0.882 


Ar3-2 


10.7 


675 


YCT8 


1174 


0.339 


0.366 


152.7 


66.6 


0.809 


Ar4 


7.4 


645 


YCT 4 


1138 


0.353 


0.377 


166.8 


147.8 


0.830 


Ar5 


7.3 


740 


YCT 9 


1174 


0.341 


0.393 


181.7 


141.4 


0.874 


Ar6 


8.8 


670 


YCT 2 


1139 


0.352 


0.381 


196.7 


150.5 


0.846 


Ar7 


8.0 


675 


YCT 11 


1128 


0.356 


0.388 


211.7 


135.0 


0.830 


Ar8 


8.2 


635 


YCT 8 


1106 


0.342 


0.361 


226.7 


150.4 


0.839 


Ar9 


7.9 


715 


YCT 3 


1154 


0.344 


0.388 


241.7 


136.6 


0.853 


ArlO 


8.1 


655 


YCT 5 


1126 


0.355 


0.381 


256.7 


149.3 


0.802 



cut was made from 15 m before to 3 h after every event that 
saturates the energy scale. 

The energy windows for event selection were our standard 
2 FWHM width (98% acceptance) and the window width in 
rise-time was set at our standard 96% acceptance. All events 
inside these windows during the entire period of counting 
were then considered as candidate 71 Ge events. Examples of 
the resultant spectra are shown in Fig. 5. It is evident that 
the energy and rise-time windows include the vast majority of 
71 Ge events. 



VII. MEASURED PRODUCTION RATE 

The number of selected events is given for each run in Ta- 
ble VII in the combination of both L and K peaks, and in Ta- 
bles VIII and IX in the L and K peaks separately. The times of 
occurrence of the candidate 71 Ge events were analyzed with 
our standard maximum-likelihood program [20] to separate 
the 71 Ge 1 1 .4-d decay events from a constant rate background. 
This is the same program that we used to analyze the runs with 
the 51 Cr source and use to analyze all solar neutrino data. The 
likelihood function is given in [1]; the only analysis changes 
from the function used in the Cr experiment were to switch 
to the 37 Ar half-life, to set the reference time to the present 
value, to slightly change the fixed solar neutrino rate to con- 
form to the current best fit for all Ga experiments of 68. 1 SNU 
[21], and to set the carryover term to zero. The latter was done 
because second extractions were made before each source ex- 
posure so there was no 71 Ge carryover from one extraction to 
the next. The number of selected events that are fit to 71 Ge is 
given in the second column of Tables VII, VIII, and IX. The 
remaining candidate events are assigned to background. The 
best fit assigns more than 200 events to 71 Ge and has a signal 
to background ratio of 1.2. For comparison, in our 51 Cr ex- 
periment, we assigned 144 events to 71 Ge and had a signal to 
background ratio of approximately 1 .0. 

Note that the first extraction Ar 1 had a strange anomaly: 
18.3 events were assigned to 71 Ge in the L peak, whereas only 
2.8 events were so assigned in the K peak. This appears to be 
just a statistical fluctuation as the counting system was func- 




FIG. 5: Upper panel: energy vs rise-time histogram of all events after 
time cuts observed in all ten exposures during the first 30 days after 
extraction. The live time is 289.0 days and 472 events are shown. 
The expected location of the 71 Ge L and K peaks based on the 55 Fe 
calibrations is shown darkened. Lower panel: the same histogram 
for the 217 events that occurred during an equal live- time interval 
beginning at day 100 after extraction. The 71 Ge has decayed away 
and is absent. The number of events outside the peaks is about the 
same in both panels as these are mainly due to background. 
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TABLE VII: Results of analysis of L- and AT-peak events. All production rates are referred to the time of the start of the first exposure. The 
combined result excludes extraction Ar 3-2. The parameter Nw 2 measures the goodness of fit of the sequence of event times [22, 23]. The 
probability was inferred from Nw 2 by simulation. 





Number of 


Number 


Number of events assigned to 


"Ge production rate 






Extraction 


candidate 


fit to 


J 'Ar source 


Solar v 


by 37 Ar source 




Probability 


name 


events 


71 Ge 


production 


production 


(atoms/day) 


Nw 2 


(percent) 


A i- 1 


9Q 


on 1 


19.4 


0.7 


1 n n+i.2 


n ha'n 
U.UOj 


OU 


Ar2 


48 


29.9 


28.7 


1.2 


10.5 + | : | 


0.048 


73 


Ar3 
Ar3-2 


69 
13 


52.9 
2.4 


51.3 

2.3 


1.6 
0.1 


14 5 +2:3 
0-8^ 


0.110 
0.273 


35 
7 


Ar4 


45 


25.4 


23.8 


1.6 


9-5-11 


0.142 


13 


Ar5 


38 


25.6 


23.8 


1.8 


11 5 +2 - 9 


0.108 


29 


Ar6 


34 


11.6 


9.7 


1.9 


6 5 +32 

U.J_ 2 .7 


0.042 


81 


Ar7 


18 


8.4 


6.7 


1.7 


6 1 +3 - 3 
u.i_ 27 

14 5 +63 


0.079 


43 


Ar8 


30 


12.8 


11.3 


1.6 


0.051 


72 


Ar9 


20 


9.3 


7.6 


1.7 


12 5 +66 


0.066 


82 


ArlO 


39 


7.4 


5.8 


1.6 


13^ 


0.151 


18 


Combined 


369 


204.1 


188.8 


15.2 


11.0^ 


0.048 


77 



TABLE VIII: Results of analysis of L-peak events. The production rate for each exposure is referred to its starting time. The production rate 
for the combined result is referred to the start of the first exposure. See caption for Table VII for further explanation. 





Number of 


Number 


Number of events assigned to 


"Ge production rate 






Extraction 


candidate 


fit to 


^'Ar source 


Solar v 


by 37 Ar source 




Probability 


name 


events 


71 Ge 


production 


production 


(atoms/day) 


Nw 2 


(percent) 


Ar 1 


20 


18.3 


18.0 


0.3 


20.2+^ 


0.094 


53 


Ar2 


23 


14.6 


14.1 


0.6 


9-3:1 


0.068 


56 


Ar3 


37 


29.0 


28.3 


0.8 


12.6^ 


0.048 


78 


Ar3-2 


10 


0.1 


0.0 


0.1 


+L1 


0.311 


8 


Ar4 


28 


16.8 


16.1 


0.8 


7 3 +2 - 5 

'"'-1.7 


0.093 


30 


Ar5 


21 


11.8 


11.0 


0.8 


4 6 +22 


0.076 


45 


Ar6 


24 


9.2 


8.3 


0.9 


3.5+ 21 


0.047 


72 


Ar7 


10 


3.1 


2.3 


0.8 


1 + ^ 


0.097 


35 


Ar8 


17 


3.1 


2.3 


0.8 


-9-5 


0.052 


74 


Ar9 


12 


5.4 


4.7 


0.8 




0.068 


51 


ArlO 


27 


0.8 


0.0 


0.8 




0.103 


48 


Combined 


219 


114.6 


107.4 


7.3 


13.2*13 


0.078 


50 



TABLE IX: Results of analysis of A"-peak events. The production rate for each exposure is referred to its starting time. The production rate for 
the combined result is referred to the start of the first exposure. See caption for Table VII for further explanation. 





Number of 


Number 


Number of events assigned to 


"Ge production rate 






Extraction 


candidate 


fit to 


J 'Ar source 


Solar v 


by 37 Ar source 




Probability 


name 


events 


71 Ge 


production 


production 


(atoms/day) 


Nw 2 


(percent) 


Ar 1 


8 


2.8 


2.5 


0.4 


2 5 +2 - 8 

-2.0 


0.040 


88 


Ar2 


25 


15.2 


14.6 


0.6 


8 9 +32 


0.092 


33 


Ar3 


32 


24.1 


23.2 


0.9 


Q 9+2.4 
—1.8 


0.100 


35 


Ar3-2 


3 


2.6 


2.5 


0.1 


i 9+O.6 
1 - i -1.0 


0.062 


79 


Ar4 


17 


8.6 


7.8 


0.8 


o o + 1.9 
-1.1 


0.075 


43 


Ar5 


17 


13.8 


12.9 


1.0 


^•'-i.i 


0.069 


64 


Ar6 


10 


3.7 


2.7 


1.0 


1 1 +L4 

1-1 -0/4 


0.081 


41 


Ar7 


8 


5.1 


4.3 


0.9 


1 7 +L4 


0.031 


95 


Ar8 


13 


9.2 


8.4 


0.8 


35^ 
-'-'-l.i 


0.116 


27 


Ar9 


8 


4.2 


3.3 


0.9 


1 3 +L3 


0.127 


22 


ArlO 


12 


5.7 


4.9 


0.8 


9-2 


0.085 


39 


Combined 


150 


91.2 


83.2 


8.0 


9 3 +1 - 3 


0.092 


34 



tioning normally and examination of the data reveals no can- 
didate events in the vicinity of the predicted K-peak position. 

As presented in [18], the counting efficiency is calculated 
based on the measured volume efficiency, the fraction of 



counting gas that is GeH4, and the counter pressure. The lat- 
ter two parameters are given in Table VI. All counters used in 
these runs have had their efficiency directly measured with ei- 
ther 71 Ge, 69 Ge, or 37 Ar. Some counters were calibrated with 
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more than one of these isotopes. 

The derived production rate of 71 Ge from the source is in 
column 6 of Tables VII, VIII, and IX. For all runs combined 
the best fit rate is H.O^-JJ atoms of 71 Ge produced per day 
by the source at the reference time. The stated uncertainty is 
purely statistical and is given with 68% confidence. Energy 
weight factors were used in this analysis in the same way as 
we analyze the solar neutrino data. Their effect on the over- 
all result is, however, quite small: without weights the rate 
decreases by 1.1%. 

Some of the systematic uncertainties that enter this result 
have been given in the foregoing. For most of the other sys- 
tematic effects we adopt the same percentage values as we 
used for the 51 Cr experiment. The results are given in Table X. 



TABLE X: Summary of the contributions to the systematic uncer- 
tainty in the measured neutrino capture rate. Unless otherwise stated, 
all uncertainties are symmetric. The total is taken to be the quadratic 
sum of the individual contributions. For comparison, the statistical 
uncertainty in the result of the 37 Ar experiment is *gjj%. 



Origin of uncertainty 



Uncertainty (%) 



Chemical extraction efficiency 






Mass of added Ge carrier 


2.1 




Amount of Ge extracted 


3.5 




Carrier carryover 


0.5 




Mass of gallium 


0.5 




Chemical extraction subtotal 


4.1 




Counting efficiency 






Calculated efficiency 






Volume efficiency 


0.5 




Peak efficiency 


2.5 




Simulations to adjust for counter filling 


1.7 




Calibration statistics 






Centroid 


0.1 




Resolution 


0.3 




Rise time cut 


0.6 




Gain variations 


+0.5 




Counting efficiency subtotal 


+3.2, 


-3.1 


Residual radon after time cuts 


-1.7 




Solar neutrino background 


0.4 




71 Ge carryover 


0.0 




Total systematic uncertainty 


+5.2, 


-5.4 



The quadratic combination of all these systematic uncer- 
tainties is ^54%. The measured production rate in the K and L 
peaks, including both statistical and systematic uncertainties, 
is thus 

p(measured) = 1 1 .O^g (stat) + 0.6 (syst) (2) 

atoms of 71 Ge produced per day. This result is not sensitive 
to any one of the extractions; e.g., if the somewhat anoma- 
lous extraction Ar 1 is deleted from the data set, the combined 
result increases by only 0.7%. 

VIII. PREDICTED PRODUCTION RATE 

For a neutrino source of activity A, it follows from the defi- 
nition of the cross section cr that the capture rate p of neutrinos 



in a material around the source can be written as the product 

p = AD(L)cr, (3) 

where D = pNofi/M is the atomic density of the target iso- 
tope (see Table XI for the values and uncertainties of the con- 
stants that enter D), and (L) is the average neutrino path length 
through the absorbing material, which in the case of a homo- 
geneous source that emits isotropically is given by 

<L>=J-f d V A \ ™L. (4) 

™'S J absorber J source 

In this last equation r SA is the distance from point S in the 
source to point A in the absorber and the source and absorber 
volumes are Vs and Va, respectively. 

The Ga-containing reactor in which the 37 Ar source was 
placed was nearly cylindrical, with a dished bottom. Based on 
accurate measurements of the reactor shape, the path length 
(L) was determined by Monte Carlo integration over the 
source and absorber volumes to be 72.6 ± 0.2 cm. The ac- 
curacy of this integration was verified by checking its predic- 
tions for geometries that could be calculated analytically and 
by noting that the measured Ga mass contained in the reac- 
tor volume agreed with that predicted by the integration. The 
sensitivity of (L) to the reactor geometry, to the position of the 
source in the Ga, and to the spatial distribution of the source 
activity were all investigated by Monte Carlo integration, and 
the uncertainty given above includes these effects. 

Based on the source activity of 409 + 2 kCi, and combining 
the uncertainty terms in quadrature, the predicted production 
rate is thus 



^(predicted) = 
atoms of 71 Ge produced per day. 



n 9 +10 

1J - y -0.4 



(5) 



IX. SUMMARY 

The measured production rate for each run is plotted in 
Fig. 6 and the predicted rate is also shown for comparison. 
The sequence of measurements fits together very well with 
^ 2 /DOF = 8.6/9 (probability of 48%), where the comparison 
is made to the combined best fit of 1 1.0 atoms/d. 

The ratio of measured to predicted production rates is 



/^(measured) 
p(predicted) 



11.0^;° (stat) +0.6 (syst) 
13 Q+ 10 



- 79+ 09 



(6) 



where the statistical and systematic uncertainties have been 
combined in quadrature. This result is nearly 2.5cr less than 
unity which has a probability of slightly more than 1%. 

To check that this unexpectedly low result is not the conse- 
quence of some experimental problem, all aspects of the ex- 
periment were carefully scrutinized, including remeauring the 
mass of Ga, verifying the position of the source within the Ga, 
and checking the extraction efficiency, the source strength, the 
counting efficiency, and the functioning of the counting sys- 
tem. No significant problems were found. 
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TABLE XI: Values and uncertainties of the terms that enter the calculation of the predicted production rate. All uncertainties are symmetric 
except for the cross section. 



Uncertainty 



Term 


Value 


Magnitude 


Percentage 


Atomic density D = pNafi /M 








Ga density p (g Ga/cm 3 ) [24] 


6.095 


0.002 


0.033 


Avogadro's number N (10 23 atoms Ga/mol) 


6.0221 


0.0 


0.0 


71 Ga isotopic abundance / ; (atoms 71 Ga/100 atoms Ga)[25] 


39.8921 


0.0062 


0.016 


Ga molecular weight M (g Ga/mol) [25] 


69.72307 


0.00013 


0.0002 


Atomic density D (10 22 atoms 7 'Ga/cm 3 ) 


2.1001 


0.0008 


0.037 


Source activity at reference time A (10 16 37 Ar decays/s) 


1.513 


0.007 


0.5 


Cross section cr [1(T 46 cm 2 /( 71 Ga atom 37 Ar decay)] [8] 


70.0 


+4.9,-2.1 


+7.0,-3.0 


Path length in Ga <L> (cm) 


72.6 


0.2 


0.28 


Predicted production rate ( 71 Ge atoms/d) 


13.9 


+ 1.0, -0.4 


+7.0, -3.0 



X. DISCUSSION 



The major purposes in making the 37 Ar source reported 
here were to develop the technology of source fabrication, to 
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FIG. 6: Upper panel: comparison of measured total production rate 
for each extraction with predicted rate. Lower panel: measured rates 
from the 37 Ar source extrapolated back to the start of the first extrac- 
tion. The combined results for events in the the L and K peaks and 
for all events are shown separately at the right and compared to the 
predicted rate. 
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FIG. 7: Results of all neutrino source experiments with Ga. The 
hashed region is the weighted average of all four experiments. 



prove that a very intense source could be made, and to elab- 
orate several techniques for source intensity measurement. 
These goals were achieved, and the source was further used 
to measure the response of the SAGE detector to 37 Ar neutri- 
nos. 

The 37 Ar source used in this experiment was made as a 
prototype for the production of a much more intense source. 
Based on the experience gained in making this source, the re- 
actor engineers for BN-600 conclude that sources in the range 
of 2.0-2.5 MCi could be made if the Ca-containing modules 
were placed in the core of the reactor, rather than in the blan- 
ket zone, as was done here. 

Since other experiments have given us great confidence in 
our knowledge of the various efficiencies in the SAGE detec- 
tor, we do not consider this experiment to be a measurement 
of the entire throughput of SAGE. Rather, we believe this ex- 
periment should be considered in combination with the other 
source experiments with Ga and interpreted as a measurement 
of the cross section for the reaction 71 Ga(v e , e~) 71 Ge. 

To this end, the results of the four neutrino source experi- 
ments with Ga given in Table I are shown graphically in Fig- 
ure 7. The weighted average value of R, the ratio of measured 
to predicted 71 Ge production rates, is 0.88 ± 0.05, more than 
two standard deviations less than unity. Although not statis- 
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tically conclusive, the combination of these experiments sug- 
gests that the predicted rates may be overestimated. 

Since 95% of the 71 Ga neutrino absorption cross section 
simply depends on the ft\/2 value for the transition from the 
ground state of 71 Ge to the ground state of 71 Ga, whose un- 
certainty is <0.5% [8], any error in the predicted rates must 
come from the contribution of the excited states. As discussed 
earlier, the Gamow-Teller strengths assigned to those tran- 
sitions were deduced from (p, n) cross sections, assuming a 
simple proportionality between (p, n) and allowed weak in- 
teraction cross sections. Yet it is known phenomenologically 
that (p, n) cross sections depend not only on the (weak in- 
teraction) Gamow-Teller amplitude, but also on a spin-tensor 
amplitude. Strong (p,ri) transitions require strong Gamow- 
Teller amplitudes, as the spin-tensor amplitude is generally a 
correction to the dominant Gamow-Teller term. In the case 
of a weak transition, however, it is possible that the spin- 
tensor amplitude dominates the (p,n) cross section. There 
are several known examples of this, e.g., the ^-forbidden Ml 
transition in 39 K— > 39 Ca [26]. In this case the Gamow-Teller 
strength contributing to beta decay is very small, yet the (p, ri) 
cross section is appreciable and attributed to the presence of 
the spin-tensor interaction. In the case of 71 Ga— » 71 Ge, the 
weak transitions to the excited states similarly could be due to 
the ^-forbidden transition amplitude of the form 1 f 5/2(1) — > 
2p3/2(p)- 



Thus, there is a theoretical uncertainty in the neutrino cap- 
ture cross section, and it is quite possible that the Gamow- 
Teller strengths to the excited states are negligible, despite the 
nonzero (p,n) cross sections [27]. As evidence for this hy- 
pothesis, we note that the weighted average of the four neu- 
trino source experiments is 0.88 + 0.05, reasonably consistent 
with R = 0.95, the value obtained if the excited state contri- 
bution were set to zero. 
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